Driven depinning of strongly disordered media and anisotropic mean-field limits 



M. Cristina Marchetti, A. Alan Middleton, Karl Saunders, and J. M. Schwarz 
Physics Department, Syracuse University, Syracuse, NY 13244 
(Dated: February 2, 2008) 

Extended systems driven through strong disorder are modeled generically using coarse-grained 
degrees of freedom that interact elastically in the directions parallel to the driving force and that 
slip along at least one of the directions transverse to the motion. A realization of such a model is a 
collection of elastic channels with transverse viscous couplings. In the infinite range limit this model 
has a tricritical point separating a region where the depinning is continuous, in the universality class 
of elastic depinning, from a region where depinning is hysteretic. Many of the collective transport 
models discussed in the literature are special cases of the generic model. 



Nonequilibrium transitions from stuck to moving 
states underlie the physics of a wide range of phenom- 
ena, from fracture and earthquake rupture to flux flow in 
type-II superconductors 0. Various dynamical models 
have been proposed in different contexts. One class of 
models, overdamped elastic media pulled by an applied 
force F, has been studied extensively. These predict a 
nonequilibrium phase transition from a pinned state to 
a sliding state at a critical value Ft of the driving force. 
Generically, this depinning transition displays critical be- 
havior as in equilibrium continuous transitions 0, with 
the medium's mean velocity v acting as the order pa- 
rameter. In overdamped elastic media, the sliding state 
is unique and no hysteresis can occur 0]. Universality 
classes have been identified, which are distinguished, for 
example, by the range of interactions or by the periodic- 
ity (or nonperiodicity) of the pinning forces. 

The elastic medium model is often inadequate to de- 
scribe many systems which exhibit plasticity (e.g., due 
to topological defects in the medium) or inertial effects 
(underdamping.) The dynamics in plastic systems can 
be inhomogeneous, with coexisting pinned and moving 
regions. The depinning transition may be discontinuous 
(first order), possibly with macroscopic hysteresis. Sev- 
eral mean-field models of driven extended systems with 
locally underdamped relaxation orphase slip have been 
proposed in the literature S H |E 13, II ■ 

In this paper, we present a model of driven disordered 
systems that encompasses several of the models discussed 
in the literature. This model incorporates the anisotropy 
of the sliding state in the plastic flow region that results 
either from flow along coupled channels oriented in the 
direction of drive (e.g., as in the moving smectic phase 
|lOl |') or in layered materials such as the high- Tc cuprate 
superconductors. We restrict ourselves to systems with a 
periodic structure along the direction of motion, such as 
charge density waves (CDWs), 2D colloidal arrays, and 
vortex lattices in type-II superconductors and consider 
only the dynamics of a scalar displacement field. This 
model uses coarse-grained degrees of freedom that are 
solid-like regions. These volumes slip relative to each 
other in the transverse dimensions, due to the presence 
of small scale defects (phase slips, dislocations, grain 
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FIG. 1: A two-dimensional realization of the anisotropic 
driven medium described in the text. Degrees of freedom 
are labeled by discrete indices £ and i, respectively transverse 
and longitudinal to the direction of the driving force, F. Each 
degree of freedom interacts with its neighbors via elastic cou- 
plings in the longitudinal direction and via viscous or similar 
slip couplings in the transverse direction. 



boundaries) at their boundaries, but remain elastically 
coupled in the longitudinal dimensions. Our current re- 
sults are for the mean-field limits, in the transverse di- 
mensions at least, with transverse viscous interactions. 
The studies carried out so far of these types of models 
for finite range interactions 0,^3 suggest that the mean 
field approximation described here may give the correct 
topology for the phase diagram, although there will be 
corrections to the critical behavior in finite dimensions. 

Anisotropic slip model. Elastic models have their lim- 
its in describing real materials. Coppersmith argued that 
elastic models with weak disorder have arbitrarily large 
strains, as large rare regions with aty pically low pinning 
result in large displacement gradients This leads 

to the breakdown of the elastic model, though the low 
density of such regions may preclude its easy observa- 
tion. Simulations, imaging and noise experiments for 
moderate and strong disorder also indicate the break- 
down of the elastic model El S [H H E ES EH. 
These calculations and observations often show a sliding 
state that consists of coherent regions moving at different 
average velocities, with extended defect structures along 
their boundaries. Furthermore, the dynamics is strongly 
anisotropic, with slip occurring preferentially where the 
shear deformations are largest, at the borders of channels 
that are on average aligned with the direction of motion. 
Guided by this work, we propose an anisotropic descrip- 
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tion of the inhomogeneous dynamics in the presence of 
defects. 

We consider a d = + d^-dimensional medium com- 
posed of c?|| -dimensional elastic channels coupled via in- 
teractions that allow for slip of the channels in the re- 
maining d± directions. The system is driven by a uni- 
form force F along one of the dy directions. A cartoon of 
a 1 -|- 1-dimensional system is shown in Fig.^ We choose 
to discretize space in both transverse and longitudinal 
directions, using integer vectors i for the dn -dimensional 
intra-layer index and £ for the d^ -dimensional layer in- 
dex. The local displacement along the direction of motion 
is (piit). Assuming overdamped dynamics, the equation 
of motion in the laboratory frame (22] is 

m) ^Kj2i^l -4^,) + F + hlYicj,} - (31) + al, , (1) 

where the dot denotes a time derivative (we have chosen 
to scale time so that the damping constant is unity) and 
(j) ranges over sites j that are nearest neighbor to i. The 
first term on the right hand side represents an elastic 
intra-channel coupling of strength K. The third term 
is the pinning force. The function Y{x) has period 1, 
and the (3^ are random phases chosen independently and 
uniformly in [0, 1). The random pinning strengths h\ are 
chosen from a probability distribution p{K). Finally, af^ 
represents an interaction (of type a) that allows for local 
slips of neighboring channels. In this paper, we assume 
a linear viscoelastic {V) stress-strain relation, 

= /* E '^™^(^ - i'WrrXt') ~ m] , (2) 

where (m) indexes the layers neighboring £. The stress 
transfer function P, 0| J^^ is generally nonlocal in space 
and/or time. We focus on a simplified version of the 
viscoelastic model where the local slip force is a purely 
viscous couphng of strength 77: ct^- = ?yX](m)(<^m ~ "^f)- 
With this coupling, Eq. is a simplified "viscous 
slip" form of the hydrodynamics of a driven viscoelastic 
medium that flows in response to large-scale shear, but 
responds elastically within the layers to long-wavelength 
compressions 0|. Gradients of the displacement trans- 
verse to the channels or layers correspond to shear de- 
formations arising from defect structures between the 
channels. Some justification for modeling the slip that 
results from moving defects as viscous couplings comes 
from work [2^ showing that that the hydrodynamics of a 
two-dimensional crystal with free dislocations is identical 
to that of a viscoelastic fluid. In contrast, gradients of 
displacement longitudinal to the drive mainly yield com- 
pressional deformations (exactly so when dy — 1.) The 
elastic response to long wavelength compressions in flu- 
ids is intimately related to the fact that compressions 
are associated with fluctuations in the conserved density. 



suggesting that the inclusion of compressional forces of 
strength K is necessary to describe the driven dynamics 
of systems with a conserved number of particles, such 
as vortex lattices. The coarse-grained model described 
by Eq. ^ allows us to investigate the competition be- 
tween elastic interaction and plastic flow in controlling 
the topology of the nonequilibrium phase diagram. 

Mean-field treatment. One mean-field approximation 
for the viscous slip model is obtained by taking the limit 
of infinite range interactions in both the transverse and 
longitudinal directions. Each discrete displacement then 
couples to all others only through the mean velocity, v = 
J2 4>i, and the mean displacement, = ^ 0^. 
We look for solutions moving with stationary velocity: 
4> — vt. Since all displacements are coupled, they can 
now be indexed by their disorder parameters (3 and h, 
rather than the spatial indices i and £. The mean-field 
dynamics is governed by the equation 

(1 + 77)0(/3, h) = K{vt -<l>)+F + rjv + hY{(t> - /?). (3) 

When if = 0, the mean field velocity is determined by 
the self-consistency condition ((jj^h))^ = v, where the 
subscript h indicates an average over the distribution of 
pinning strengths p{h). When K 0, the mean field 
velocity is found by imposing {(j>{h, f3) — vt)h.p — 0. 

It is useful to review the case where 77 = and if 7^ 0. 
In this limit, Eq. Q reduces to the mean field theory 
of a driven elastic medium worked out by Fisher and 
collaborators [25j |. For the piecewise harmonic pinning 
y(0) = 1/2 — for < (/) < 1, no moving solution ex- 
ists for F < Ft = { 2(K+h) "^^^^ above threshold the 
mean velocity has a universal dependence on the driv- 
ing force, with v ^ {F — Ft)^ ■ The critical exponent 
f3 depends on the shape of the pinning force: /? = 1 for 
the discontinuous piecewise harmonic force and /3 = 3/2 
for generic smooth forces. Using a functional RG expan- 
sion in 4 — e dimensions, Narayan and Fisher |25 j| showed 
that the discontinuous force captures a crucial intrinsic 
discontinuity of the large scale, low-frequency dynamics, 
giving the general result /? ~ 1— e/6-|-0(e^), in reasonable 
agreement with numerical simulations in two and three 
dimensions H^]. For simphcity and to reflect the 
"jerkiness" of the motion in flnite-dimensional systems 
at low velocities, we use piecewise harmonic pinning. 

When J] > 0, the nature of the depinning differs qual- 
itatively from the r] = case in that hysteresis in the 
dynamics can take place. This can be shown starting 
from the 77 = solutions w^=o(^)i with Vr/^o a well- 
deflned function of F. The solution for general 77 can be 
found by substituting the effective drive force G — F-\-riv 
for F in the Vr]=o{F) relation and scaling velocity down 
by 1 -|- 77. The linear transformation F — G — rjv then 
gives the general v{F) curve. One result is that when 
1 + 77^^ < vaax.F{dvri=o / dF) , the v vs. F curve is multi- 
valued. Two cases can now be considered. The viscous- 
elastic case K ^ Q has quite different behavior from 
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FIG. 2: Velocity curves and phase diagrams for the mean- 
field slip model. Part (a) shows the velocity ii vs. force F 
curves for the purely viscous case for a broad distribution 
of pinning, p(K) = e~*^ , for two values of viscous coupling rf. 
Part (b) shows the phase diagram from (a), with curves curves 
indicating jumps in the adiabatic response (i.e., dv/dF — > 
oo) when the force F is increased from F = {F^ (rj)) and 
decreased from a sliding state {F^{ri)). The critical point 
{r]c, Fc) « (9.61, 0.278). Part (c) shows example v vs. F curves 
when there is elasticity, K = 1, for piecewise linear pinning 
forces and p{h) = e"**. The curves are again continuous for 
ri < rfc, but have a depinning transition at Ft, with critical 
behavior identical to the rj — case. When r) > rjc, the 
depinning is abrupt and the v-F curves are hysteretic. Part 
(d) summarizes the if = 1, p{h) — e"^ results in a phase 
diagram, indicating depinning from an initial pinned state at 
F^ = Ft and repinning from a sliding state at a lower force 
F^{ri). The point (r^c, Ft) ^ (3.74, 0.298) is a tricritical point. 



the purely viscous case if = 0, due to the distinct self- 
consistency conditions imposed by elasticity (or conser- 
vation in the flow direction.) 

For the purely viscous case K = Q (the r = limit of 
the model in Ref.0), periodic solutions of period T{h, G) 
are found for a single particle moving in a pinning force of 
strength h under the effective drive G. Self-consistency 
implies that v = {\T{h,G)]~^)h, giving 
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The behavior now depends on the shape of the disorder 
distribution p{h). For narrow distributions, that vanish 
below a value /imin > 0, no moving solutions exist for F < 
^inin/2. In the extreme case of uniform pinning, p(h) = 
d{h—ho), the depinning is hysteretic for all nonzero values 
of T], with a jump from w = to finite v at F = ho/2. For 
the more general broad distribution with nonvanishing 
support at /i = 0, no pinned solution exists for F > 0: 
while the bulk of the degrees of freedom are pinned at 
small F, weakly pinned (/) can respond to small drives. 
Typical v vs. F curves for p{h) = are shown in Fig. 
2(a). The lines in Fig. 2(b) represent the critical forces 



where there is a macroscopic jump in the average velocity, 
indicating dynamic hysteresis. For rj > rjc and F-^ > 
F > Fi "fast" and "slow" sliding states can coexist. The 
critical point (r]^, Fc) in Fig. 2(b) is in the universality 
class of the liquid-gas transition (and of the field-driven 
random field Ising model 0| .) Near the critical point the 
mean velocity has universal scaling, v — Vc ^ {'H^'Hc)^^^ ^ 
{F-Fcf, witht;^ = v{Fc,rjc), and Pmf = 1/2, Smf = 3. 
The point (77c, Fc) survives in finite dimensions 11], with 
exponents distinct from their mean field values. 

For finite long-time elasticity, i.e., when K ^ 0, the 
behavior changes dramatically. The elastic forces or par- 
ticle conservation enforce a uniform time-averaged veloc- 
ity for all degrees of freedom. In this case the long-time 
uniform-w solution to Eq. © is (j){h, /3, t) = vt + <j)[h, /3, t) 
with <j>{h,P,t) of period 1/v in t. The self-consistency 
condition is {4>{f3,t)) — 0. The explicit solution can 
be obtained for the piecewise linear pinning force. (In 
the mean-field limit the nature of the phase diagram 
depends some on the shape of the pinning potential.) 
There is now a pinned phase that is stable |2^ for 
F < Ft = i ^rr'^ , )h- One then solves 
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(K + h) e(^+M/[(i+';)^'] -l/h 



(5) 



to find the mean velocity in the sliding state (0 < 6F = 
F - Ft), where {ijc + 1)"^ = i (K+hr^ )h- The velocity- 
force curves and a phase diagram are shown in Fig. 2(c,d) 
for p{h) = e^^ . There is now a tricritical point at 
{rjc, Fc = Ft). For rj < rjc, a continuous depinning tran- 
sition at Ft separates a pinned state from a sliding state 
with unique velocity w ~ (1 -f rjc) , giving [3=1 
in MFT. In finite dimensions, this transition is likely to 
remain in the same universality class as the depinning 
of an elastic medium (77 = 0): numerical studies and 
analysis by Schwarz and Fisher |l2l | of a model with lo- 
cal viscous-like terms show that there is no hysteresis for 
small slip coupling and that the depinning transition ex- 
ponents are the same as without slip coupling. In our 
mean- field example, the linear response diverges at 7]c, 
v{rj = rjc) ~ 1/ ln(F — Ft). For rj > rjc there is hysteresis 
with coexistence of stuck and sliding states. 

Relations to other models. In the absence of viscous 
coupUng {rj = 0), Eq. ^ reduces to the conventional d\\- 
dimensional phase-only model of driven CDWs. When 
K = Q and rj is finite, the model describes shear defor- 
mations in a viscous fluid in a disordered background Q ■ 

The mean-field limit of the anisotropic slip model for 
finite K maps onto a periodic-pinning version of the 
infinite-range version of the stress overshoot model stud- 
ied by Schwarz and Fisher [^^3- The dynamics of crack 
fronts in brittle materials can be dominated by local in- 
ertial effects, in which the motion of a crack segment 
creates a transient stress on other segments, the stress 
overshoot, resulting in a net forward motion beyond the 
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local minimum of the pinning. In the infinite-range limit 
the stress transfer yields a global coupling of strength M 
between the mean displacements at two different times, 
analogue to the effective force per unit velocity 77 in our 
model. The behavior we conjecture for finite dimensions 
is identical to that obtained by Schwarz and Fisher. 

The viscous model (77 ^ 0) with finite K is also closely 
related to a model of sliding CDWs that incorporates 
the coupling of the CDW to normal carriers by adding 
a global velocity coupling 0, to the conventional 
Fukuyama-Lee-Rice model. The precise relationship is 
obtained by considering a mean-field limit — > 00 with 
fixed d\\. In this limit, each degree of freedom is coupled 
equally by slip to all other layers, but elastically coupled 
only to its nearest neighbors in the longitudinal direction. 
The equation of motion is then 
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>^)+F + r;« + /i^y(</>^ -/?;). (6) 



Eq. © is identical to the equation of motion of a purely 
elastic CDW with friction 1 + -q and an effective driving 
force F+rjv. One can then obtain the velocity- force char- 
acteristics consistent with Ref. '9' simply by translating 
the analytic results |2,5] for CDWs. Near threshold, we 
obtains - (F-t-r/w-Fr)'^, with /3 = l-(4-d||)/6+0[(4- 
The transition is always hysteretic: the medium 
depins at F"^ = Ft when the force is ramped up and re- 
pins at a lower value, F^ fa Ft — (^^ - 1^ (1^) ' 
when the force is ramped down from the sliding state. 
This hysteresis for all 77 > appears to be a consequence 
of the global nature of the coupling, in contrast with the 
results of Schwarz and Fisher for local slip coupling ■ 

One can include phase-slip couplings a^^ periodic in 0. 
The anisotropic model can describe the driven dynamics 
of intrinsically layered systems, such as vortex lattices in 
cuprate materials. Hysteresis depends on the strength of 
the correlations in the (weak-coupling) layered direction, 
as seen in simulations 17] . An anisotropic CDW model 
was studied in this context by Nattermann and Vinokur 
0. These authors modeled the transverse slip as a non- 
linear coupling that is periodic in the local CDW phase 
differences; this allows for phase slips, i.e., the collapse 
of the CDW amplitude between coherent regions ^ . 

Summary. We have proposed that the plastic flow 
of extended systems in a disordered background can be 
modeled using a general anisotropic approach that in- 
cludes many specific physical systems as special cases. 
We have discussed the behavior of this model in a mean- 
field limit. When the degrees of freedom are allowed to 
move at distinct velocities, a pinned phase is generally 
not present, though there is hysteresis in the current- 
drive relation. Even when fluid-like shear takes place, 
particle conservation gives a sharp depinning transition 
in flow that takes place along channels. One of our future 
goals is to establish the connection between the param- 
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